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DISPLACEMENT EFFECT OF A THREE-DIMENSIONAL BOUNDARY LAYER 1 

By Fbanklin K. Moobb 


SUMMARY 

A method is described for determining the “displacement 
surface ” of a known three-dimensional compressible boundary- 
layer flow in terms of the mass-flow defects associated with the 
profiles of the two velocity components parallel to the surface. 
The result is a generalization of the plane flow concept of dis- 
placement thickness introduced in order to describe how a thin 
boundary layer distorts the outer nonviscous flow. 

The height of the displacement surface above the body surface 
for flow about a yawed infinite cylinder is shown to be equal to 
the height characterizing the mass-flow defect of the chordwise 
velocity profile. The displacement-surface height is shown to 
differ, in general, from that associated with the resultant mass- 
; flow defect, even at stagnation points of the secondary flow. 
Numerical values are found for the known three-dimensional 
boundary-layer flow about a cone at a small angle of attack to 
a supersonic stream. 


INTRODUCTION 

The boundary layer established in the flow of a slightly 
viscous fluid about a body is normally considered an isolated 
region wherein the effects of viscosity predominate and 
outside of which the motion of the fluid is governed by the 
laws of nonviscous motion. For large Reynolds numbers, 
the boundary layer is assumed to be so thin that the non- 
viscous portion of the flow occurs as though there were no 
boundary layer. This assumption is strictly correct in the 
limit of infinite Reynolds number. For large but finite 
Reynolds numbers, the growth of the boundary layer causes 
the stream to be deflected away from the body surface. 

This displacement effect of the boundary layer on the 
nonviscous flow may properly be determined from the behav- 
ior of the boundary layer itself, as established either by 
experiment or by solution of the Prandtl boundary-layer 
i equations for laminar flow. 

It does not follow, however, that this revised outer flow 
may properly be used in conjunction with the Prandtl 
equations to yield an improvement in the boundary-layer 
calculation. Such an improvement may be obtained only 
by use of a new set of equations that take into account the 
variation of pressure across the boundary layer. This varia- 
tion is neglected in the Prandtl equations. (See Alden’s 
iterative solution for incompressible flat-plate flow, ref. 1.) 


The customary definition of displacement thickness (ref. 2) 
is applicable to two-dimensional flow and is expressed in the 
following equation: 


j: 


pu dy=piUi(h—b*) 


( 1 ) 


where h is some location well outside the boundary layer 
(fig. 1) at which pu—piu, and beyond which the flow may be 
considered nonviscous. (A list of symbols is provided in the 
appendix.) Under the assumption of an extremely thin 
boundary layer, h is so small that p,u, may be taken as the 
evaluation at the body surface (y=0) of the nonviscous flow 
obtained by neglecting the presence of the boundary layer'. 
Equation (1) equates the actual mass flow near the surface 
with the mass flow which would be associated with a non- 
viscous flow that terminates at 5* rather than extending to 
the wall. ' Thus, the nonviscous portion of the flow behaves 
as if it occurred in the presence of a solid boundary given 
by the displacement thickness 8*(x). Equation (1) may be 
solved for 8*: 

'-LX'Sd* 


Ordinarily, theoretical boundary-layer solutions for pu/piUi 
asymptotically approach 1 for large values of Reynolds 
number based on y. Therefore, displacement thickness is 
often defined as follows: 



i Supersedes NAOA TN 2722, "Displacement Effect of a Three-DImenslonfll Boundary Layer," by Franklin K. Moore, 1952. 


367 



368 


REPORT 1124 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Figure 2. — Three-dimensional boundary-layer flow (Cartesian 
coordinates). 


For three-dimensional boundary-layer flows, two lengths 
characterizing mass-flow defects may he defined in terms of 
the profiles .of the two velocity components tangential to the 
surface (fig. 2), 

(2) 

and it is not dear which, if either, defines a “displacement 
surface” that properly describes the extent to which the non- 
viscous flow is deflected by the boundary layer. 

Of course, it is expected that there does exist a displace- 
ment surface for such flows. The analysis that follows 
shows that such a surface may he described by a defining 
equation more fundamental than equation (1) . 

The velocity and density profiles are assumed to be known 
for the three-dimensional flows under consideration. Cases 
for which this is true include: the laminar boundary layer 
on yawed infinite cylinders, treated by Prandtl (ref. 3), 
B. T. Jones (ref. 4), and Sears (ref. 5); and the laminar 
boundary layer on a cone at a small angle of attack to a 
supersonic stream (ref. 6). The displacement effect of these 
flows will he treated specifically. 

The investigation was conducted at the NACA Lewis 
laboratory in February 1952. 

THEORY 

DEFINING RELATION FOE DISPLACEMENT SURFACE 

The boundary-layer solution (assumed known) yields a 
certain distribution of velocity Vt,i{x,z) normal to the body 
surface at the outer edge h(x,z ) of the boundary layer, where 
p/u, and w may be taken essentially equal to pi,Ui, and Wi 
(see fig. 2) . Under the assumption that the nonviscous flow 
is altered only slightly by displacement, the most direct 
way to compute this effect would be to suppose that the 
nonviscous equations hold for y>h(x,z) and to impose the 
following boundary condition on the normal velocity v of in 
the outer flow:. • 


However, since the boundary conditions usually encountered 
in nonviscous flow specify an impermeable surface, it is 
convenient to recast the boundary condition given in equa- 
tion (3) in answer to the question: What impermeable 
surface y=A(x,z ) would deflect a nonviscous fluid in such a 
way as to produce a normal velocity satisfying condition (3) ? 
This fictitious surface may he called the displacement surface. 

Because A(x, z) is imagined to be a stream surface, it is 
necessary to specify that at y= A the resultant velocity 
vector ( u , v, w) be tangent to the surface y— A. Thus, at 
y= A the ratio of v to the magnitude of q, which is defined as 
the vector (u, w), must be set equal to the slope of the surface 
y=A, measured in the direction of the vector q, or, equiva- 
lently, equal to the component in the direction of q of the 
vector grad A. In vector notation, therefore, the normal 
velocity v which would be produced in a nonviscous fluid at 
an impermeable surface y= A is q • grad A. The vector q 
may be obtained by evaluating the velocity vector of the 
unrevised nonviscous flow at the body surface (y— 0), under 
the related assumptions that the velocity vector varies only 
slightly over distances of the order of the actual boundary- 
layer. thickness and that the revision required to take account 
of displacement is slight. Thus, at y=A (see fig. 3), 

v 0 f = £i • grad A 

• 

The increment in v of between A and h is approximately 
(h-A)bD 0 ,[by; again, a thin boundary layer is assumed and 
only the first term in a Taylor’s series is used. 

To the order of approximation contemplated in this 
analysis, bv of /by .m&y be obtained from the unrevised non- 
viscous flow evaluated at y= 0. Thus, the fictitious im- 
permeable surface A would produce, at y=h, a normal 
velocity (see Tig. 3) 

o 0 / = 2 l .gradA-Kfc— A)(^_ fl 

Introducing this result into the boundary condition given in 
equation (3) yields the defining relation for A(x, z) 

(»«)»-*= (h~ A) o + gi • grad A ' (4) 

The boundary-layer solution yields 

EXPRESSION OF A IN TERMS OF MASS-FLOW DEFECTS 

The displacement surface A may be related to the mass- 
flow defects (eqs. (2)), which characterize the boundary 
layer, as follows: In a Cartesian coordinate system (fig. 2), 
the equation of continuity for both,the boundary-layer and 
nonviscous flow is 

bpv bpu bpw ... 

by ~ bx bz W 

Under the Prandtl boundary-layer assumptions, the Cartesian 
equations of motion, and hence equation (5), may bo applied 
in an orthogonal curvilinear coordinate system in which the 


v 0 f=v b i{x,z) at y=h{x,z) 


( 3 ) 
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surface of the body is given by y=0, provided tbe radius of 
curvature of tbe body is everywhere large as compared with 
the boundary-layer thickness. Integrating equation (5) 
across the boundary layer yields 


(p Vbl) tr-A— 




= /o [cir P u )+^ (piWi— pio)J dy— 


IX 


. bpiWi 


2>z 


2)2 


) dy 


or, inasmuch as h is outside the boxmdary layer where 
pw=Pi«i and pw—piVJi, 


Pi (,Vbi)p-h—~h (piUiS z )+~ (piwrf.) ( 6 ) 


condition must be imposed at the location A =S X , at least 
near the stagnation point. 



Figure 3. — Displacement surface in three-dimensional flow (Cartesian 

coordinates) . 


where S x and S x are the mass-flow defects defined in equations 
(2). For the nonviscous flow, v o/ =0 at the body surface 
(y=0), and equation (5) becomes 


Pi 


r\ /2>Pl«l, 

: 2>Pl10l\ 

A-o V dr 1 

h dz ) 


( 7 ) 


Introducing equations (6) and (7) into equation (4) yields 
[piiti(A— +jjj^ [piWi(A— 5 t )] = 0 (8a) 


In a Cartesian coordinate system, the displacement surface 
is related to the mass-flow defects by equation (8a). When 
cases arise for which other coordinate systems must be used, 
the following generalization in vector form of equation (8a) 
may be used: 


div J-pi^A — Jo (pigi— pq)dyj=0 (8b) 

where y is the distance normal to the body surface and where 
the divergence operator involves differentiation only with 
respect to the two coordinates parallel to the body surface. 

EXAMPLES 

Plane flow. — Equation (8a) may be integrated to yield, 
for plane flow (d/bz=0), 

A =S*+t|- (9) 

piUi 

where K is a constant of integration. 

The appearance of this constant means, in general, that 
the revised boundary condition on the nonviscous flow near 
the wall may be applied at any surface in the boundary-layer 
region; for example, along the wall itself. Of comae, if 
there is a stagnation point on- the body where vanishes, 
then K must be taken equal to zero and the revised boundary 


Stagnation point of secondary flow. — In the vertical plane 
of s ymm etry of the flow about a body of revolution at an 
angle of attack, the circumferential velocity component 
vanishes in the boundary layer as tvell as in the nonviscous 
outer region. Such stagnation of the secondary motion 
would occur in a variety of cases, in particular, wherever an 
essentially three-dimensional flow has a plane of symmetry. 
When 101 is designated as the component of secondary flow 
(circumferential velocity for a body of revolution) at the 
outer edge of the boundary layer, and equation (8a) is 
written in the form 

^^piiii(A— 5i)J— (-101 q^|^pi(A — 5,)J+Pi(A— 5.) 0 (10) 

it is clear that A cannot equal in general, even if both Wi 
and the boimdaiy-layer profile of 10 vanish, unless bw 1 fdz 
also vanishes. 

Flow about yawed infinite cylinder. — In the flow about a 
yawed infinite cylinder, there is a spanwise boundary-layer 
velocity profile and an associated spanwise mass-flow defect. 
(See refs. 3 to 5.) If x is taken in the chordwise direction 
(fig. 4), the entire flow depends only on x. Thus, derivatives 
with respect to the spanwise coordinate 2 vanish; and, hence, 
from equation (8a), the plane-flow result (eq. (9)) applies. 
Accordingly, the spanwise mass-flow defect represented by 
5 X does not enter into the determination of the displacement 
surface. 



Figure 4. — Coordinate system for yawed infinite cylinder. 
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FiatrRE 5. — Coordinate system for cone at angle of attack. 


Supersonic flow about cone at small angle of attack. — 
When the coordinate system shown in figure 5 is used, 
equation (8b) becomes 

sin 8 ^ (A — £pi 10 i(A — 5?)J=0 (11) 

where 



Because the outer flow is conical, pi, u u and Wi are functions 
only of <p ; and equation (11) may be written 

prill sin 6 ^ jjc (A— j^piWi (A— 5 P )J= 0 (1 3) 

In the case of conical outer flow, the associated boundary- 
layer profiles show similarity of the Blasius type in meridional 
planes (see ref. 7 or 8). Thus, in a meridional plane, A, 8 X , 
and 5 are proportional to -\fx. Incorporation of this infor- 
mation into equation (13) gives 

^ piU-i sin 0(A— 5,)+^ j^PiWi(A— 5 P )J=0 (14) 

For a cone at small angle of attack a, u x is nearly equal to 
u, the velocity on the cone surface at zero angle of attack. 
The quantities pi, A, 5„ and vary only slightly with angle 
of attack, whereas 

Wi=cmA 3 sin <p (15) 

where A 3 depends only on the cone vertex angle and the flow 
Mach number and is defined in reference 6 as follows: 

, z 2x 

— — — — = : a 

U U 81IL 6 

The quantities z, x, and u are in the notation of reference 9, 
wherein they are tabulated. 


To the first order in a, the substitution of equation (16) 
into equation (14) yields the result 

A=5j-{-— . „ Axis?— 8^) cos tp (16, 

o SUL u . 

The analysis of reference 6 yields the values of S x and 8 V . 

Clearly; A differs from 8 X in the plane of symmetry <p= 0, ir, 
where the circumferential velocity w vanishes. It might, 



Figube 6. — Proportional rate of increase of displacement thickness 
with angle of attack on cone. 
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however, be noted that in the plane A=5 I . Figure 

6 shows the proportional rate of increase of displacement 
thickness A with angle of attack in the plane of symmetry 
<p—ir at zero angle of attack. The corresponding rate of 
increase of the mass-flow defect S x is shown for comparison. 
These curves are obtained from equation (16) and the 
results of reference 6. The sine of the semivertex angle is 
introduced as a factor primarily to permit presentation of 
the case 6 = 0 as the limit of an indeterminate form. For a 
stream Mach number of 2, the change in displacement thick- 
ness appears to be of the order of 50 percent larger than the 
change in the meridional mass-flow defect. 


CONCLUDING REMARKS 

The foregoing analysis deals only with the displacement 
effect of a known boundary layer on the nonviscous outer 
flow, and hence can be applied only if the boundary-layer 
behavior has been determined either theoretically or experi- 
mentally. The latter approach might possibly find appli- 
cation in the correction of nozzle contours for boundary-layer 
development. In this connection, perhaps, it should be 
noted that the analysis is not restricted to laminar flows. 

Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 6, 195 2 


APPENDIX— SYMBOLS 


The following symbols are used in this report: 

Aq function of cone angle and Mach number 

(eq- (15)) 

h height above body surface at which p, u, 

w~pi, u u w h and beyond which nonviscous 
equations apply (eq. (3)) 

M 0 stream Mach number 

2 velocity vector composed of components paral- 

lel to body surface u, w 

u velocity component in r-direction 

u meridional velocity component at surface of 

cone at zero angle of attack 
v velocity component in y-direction 

v b i boundary-layer solution for velocity normal to 

surface, evaluated at outer edge of boundary 
layer 

v of nonviscous solution for velocity normal to 

surface, evaluated near the surface 
w velocity component in z-direction 


:} 

coordinates in body surface 


V 

coordinate normal to surface 


a 

angle of attack 


A 

.height above body of displacement surface 


5* 

displacement thickness in plane flow 


8 : 

length characterizing mass-flow defect 

of 


■u-profile (eqs. (2) and (12)) 


8 , 

length characterizing mass-flow defect 

of 


10 -profile (eq. (12)) 


8 r 

length characterizing mass-flow defect 

of 


10 -profile (eq. (12)) 


e 

semivertex angle of cone 


p 

density 


<p 

angular coordinate (fig. 5) 



Subscript 1 denotes evaluation of nonviscous flow at body 
surface, taken equivalent to conditions at outer edge of 
boundary layer of infinitesimal thickness. 
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